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Abstract 

In the framework of Generalised Parton Distributions, we study the hehcity- 
dependent and independent cross sections measured in Hall A and the beam spin 
asymmetries measured in Hall B at Jefferson Laboratory. We perform a global fit 
of these data and fits on each kinematical bin. We extract the real and imaginary 
parts of the Compton Form Factor 7i under the main hypothesis of dominance of 
the Generalised Parton Distribution H and twist 2 accuracy. We discuss our results 
and compare to previous extractions as well as to the VGG model. We pay extra 
attention to the estimation of errors on the extraction of TC. 

Introduction 

Since it was realised that Generalised Parton Distributions (GPD) were reachable through 
an harmonic analysis of the Deeply Virtual Compton Scattering (DVCS) process [1], [2], 
the study of the connection between GPDs and DVCS has been a very active field of 
research, concerning theoretical developments (see the reviews |3], [5] and [6]) as well 
as experimental ones ([7],[8],[9],[T0],[IT],[l2],[l3],[Tl]andfT5]). After the first dedicated 
experiments and at the beginning of the experimental GPD program, it is already worth 
trying to extract GPDs from measurements. 

The present work addresses this question, and illustrates it with recent JLab data, 
namely beam spin asymmetries (BSA) [13] and helicity-dependent and independent cross 
sections [12]. These data offer a large kinematic coverage and a fine kinematic binning, 
which are interesting features for our purpose. However, the methods we use and the 
conclusions we come to are presumably not restricted to JLab kinematics, and may be 
of interest for other experimental set-up. Furthermore, this is one of the first global fits 
of these measurements, and such fits are necessary to the completion of the experimental 
GPD program. 

*herve. moutarde@cea.fr 
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The harmonic analysis of ep ep'y cross-sections has so far rehed on the 2002 work of 
A.V. Behtsky, D. Miiller and A. Kirchner [16]. In this formahsm, the interference between 
the Bethe-Heitler (BH) and DVCS processes was treated with a leading order approxima- 
tion of the BH part. This assumption was removed by A.V. Belitsky and D. Miiller in [IT] 
in the case of a spinless target, and by P.A.M. Guichon and M. Vanderhaegen in the case 
of a proton target [18]. In all the following, we will use the expressions from the latter. 

The first section of this paper describes JLab data, the difference of the evaluation of 
the ep — >• ep7 cross sections by [TH] and [IH], and outlines our hypothesis. In the second 
part we explain our fitting strategy, which is twofold : we perform a global fit of the selected 
JLab data and fits on each kinematic bin (which we will be refering to as local fits). In 
the third section we systematically compare the results of global and local fits, and keep 
the output of the global fit as our best solution. We then discuss the extracted values of 
7i and compare them to previous studies. 

1 Preliminary analysis 

In this study, is the standard Bjorken variable, the virtuality of the initial photon 
and t the square momentum transfer between initial and final protons. 

1.1 Description of the selected JLab data 

C. Munoz-Camacho et al. [12] published hehcity-dependent and heli city- independent cross 
sections for = 0.36 and t between -0.33 and -0.17 GeV^. Helicity-independent cross 
sections are given at = 2.3 GeV^ only, while helicity-dependent cross sections are 
measured for between 1.5 and 2.3 GeV^. Data come in 12 bins (3 values of and 4 
values of t). Helicity-dependent cross sections are given with relatively large uncertainties 
(from 20 % to 100 %, depending on the bins and on the value of the angle between leptonic 
and hadronic planes) while helicity-independent cross sections are usually accurate at the 
5 % level. 

F.-X. Girod et al. [H] released BSAs over a wide kinematic range [xb from 0.11 to 0.58, 
from 1 to 4.8 GeV^ and t from -1.8 to -0.09 GeV^, described by 62 mult i- dimensional 
bins (5 values of t, 5 values of xb and 4 values of Q^). Their accuracy is 25 % on average, 
ranging from 5 to 100 % depending on the bins and on the value of the angle between 
leptonic and hadronic planes. 

These two sets of measurements have kinematic configurations in common, allowing us 
to perform cross-checks by evaluating a BSA through the ratio of helicity-dependent to 
helicity-independent cross sections. This simple exercise shows no unexpected discrepancy 
between Halls A and B measurements. This implies that both sets of data are consistent 
and can be used in a global fit. 
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1.2 DVCS at leading twist 

Four GPDs E, H and E appear at twist 2, but the cross sections depend on the 
Compton Form Factors (CFF). The convention of is used to define the CFFs : 

T= r dxF{x,U)(- — - . ^ . ] {F = HoiE) (1) 

T = f dxF{x,^,t)(- ^ - + - ^ A {F = HotE) (2) 
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and ^ = xb ^^^w^' generahsed Bjorken variable [IE], [18]. The complex integration 

kernel yields a real and an imaginary part to the CFFs : 

ReJ^ = 'P j dxF{x,^,t) ( {F = HotE) (3) 



and 



7r(F(e,e,t)-i^(-e,e,t)) {F = HoiE) (4) 



ReT = 'P j (^~^ " {F = HotE) (5) 

W = n(^F{^,U) + F{-^,U)) {F = HotE) (6) 

where the symbol V denotes the principal value of the integral. 

The reak and imaginary parts of a CFF are related by dispersion relations due to 
analycity properties (see [19], [20] or [2T]). However the unknown substraction (the D- 
term [22]) and the limited kinematic range of our data make this constraint rather weak. 
In this work we consider the real and imaginary parts as independent (local fits) except 
in the case where we use an explicit parametrisation of the GPD H to compute the CFF 
(global fit). 

1.3 The formalism of P. A.M. Guichon and M. Vanderhaeghen 
and its consequences 

The recent work of P.A.M. Guichon and M. Vanderhaeghen (GV) gives analytic expression 
for the ep — > epy cross section ; these formulae are embodied in a Mathematica package 
|18j . which we used to build a C++/ROOT library. 

The ep — > epj cross section is classicaly divided into three parts, namely BH, VCS 
and interference between the BH and VCS cross sections. The contribution of the BH 
amplitude to the interference is treated exactly. This introduces some differences with 
respect to the well-known BMK expressions and two new important qualitative features to 
the discussion of the measurements under scrutiny. Let us mention before an important 
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point for actual computation of cross sections : even if the GV expressions rely on the 
same tensorial decomposition of the hadronic VCS tensor than BMK expressions, the GV 
kinematic conventions follow those used in early versions in the VGG code detailed in [23] . 
The BMK and GV formalisms thus use different definitions for the angle between the 
hadronic and leptonic planes, and these definitions do not follow the Trento convention 
[21] chosen to describe the measurements : ^Trento = vr — 0bmk and ^Trento = 27r — 0. 

Concerning Hall A measurements, the BMK formalism restricted to twist 2 asserts that 
helicity-dependent cross sections write : 



d^^ d^^ 



C.s,„^ /m(H+^(l + f)H-^f£) (7) 



XB Fi' AM^Fi 



where d'^^ = dQ'^dxBdtd(j), Fi and F2 are the Dirac and Pauli form factor, M the proton 
mass and Ci is a constant irrelevant for our purpose. Since the Q^-dependence of this cross 
section is factorised, this expression allows a study of scaling without having to disentangle 
the different CFFs as in [12]. 

With the exact equations of [18], an helicity-dependent cross section does not have such 
a simple form. It writes : 



C2 sin (f) Ira (h + cs £ + H + + 



where C2 is a constant irrelevant for our purpose. The dots stand for power-suppressed 
contributions. The test of scaling is more involved now since the coefficients cs, Cf^ and 

do depend on Q"^. At given xb and t, the coefficients of Eq. ([s]) ce, Cf^ and vary 
respectively by 20 %, 6 % and 38 %. If we only fit the combination of CFFs appearing in 
Im{. . .) in Eq. ([s]), the kinematic Q^-dependence of C£, and may appear as a scaling 
deviation of the same magnitude. 

Differences also arise in the expression of a BSA. Its dependence on the angle between 
the leptonic and hadronic planes takes the following form : 

BSA = asin0 + 6sin20 

1 + c cos (f) + d cos 20 + e cos 30 

where a = 0{Q-^), b = 0{Q-^), c = 0{Q-^) d = 0{Q~'^) and e = 0{Q-^) are real 
numbers. In the BMK picture, the coefficients b, d and e are higher-twist contributions. 
As a straightforward consequence, we see that the 90° asymmetry is no longer proportional 
to the imaginary part of a linear combination of CFFs. 

Moreover the coefficient c of Eq. ^ now depends on the imaginary part of CFFs, and 
not only on the real part as in [16]. This precludes a clean separation of the real and 
imaginary parts of 

through the distinct measurements of the helicity-dependent and helicity-independent cross 
sections as was done previously [12] . 
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1.4 Main assumptions 



P.A.M. Guichon and M. Vanderhaeghen provide exact analytical expressions of the ep 
ep7 observables as functions of the DVCS amplitudes. The latter can be written in terms 
of GPDs through the usual twist expansion [1]. In this work we restrict ourselves to the 
twist 2 approximation. This is a reasonable assumption, since C. Muiioz-Camacho et al. 
[12] claimed the observation of early Q^-scaling. 

As we are considering experiments on a proton target, we neglect H and E {H- 
dominance) for the following reasons. Firstly, this is supported by kinematics : in Eq. (|8]), 
the coefficient C£ varies between 0.05 and 0.28, Cf^ between 0.3 and 1., and between 
2. 10~^ and 0.015 for the 52 kinematic configurations of Hall B data satisfying < ^• 

Secondly, for small t and ^, we expect ^ to be close to ^ i.e. |. Thirdly, we can check, 
for instance thanks to the VGG model [23], [25], [26] and [27], that the relative sizes of 



ImS and ImTi to ImT-C are similar : varies between 0.21 and 0.92, and between 
0.13 and 0.91 for the same set of 52 Hall B kinematic configurations. This indicates that 
the hierarchy between the kinematic coefficients refiects the hierarchy of contributions to 
the interference. 

Thus assuming i7-dominance, we may hope to extract information on H from BSAs or 
helicity-dependent cross sections with a systematic error of 20 % to 50 this approxima- 
tion being better at small t. The advantage of this approach is the dramatic decrease of the 
number of degrees of freedom involved in fits. M. Guidal indeed showed in a recent work 
[28] in the same kinematic region that it is not possible to extract sensible information 
about the real and imaginary parts of H, £, 7i and £ by direct fits of hehcity-dependent 
and independent cross sections. More specifically, keeping only the dominant coefficients 
in Eq. ^ gives the minimal functional form : 

, a sin 6 , , 

BSA ~ — 10 

1 + c cos 

A direct fit of BSAs to this reduced expression on each (xB,Q^,t)-bin, along the lines of [H], 
shows that the coefficient c is compatible with (while with a marked trend to negative 
values) within error bars for 25 bins over the 52 bins for which |4 is less than \. Extracting 



the real part of the CFFs contained in the coefficient c of Eq. (10) is thus a demanding 
task. 



2 Fitting strategies 

The possibility to study GPDs in DVCS rests on factorisation theorems [29], which require 
a small value of In the following, we restrict ourselves to kinematic configurations for 

which ^ <\- 

direct test of the ff-dominance assumption with the VGG model gives an upper bound of 25 %. 
This is comparable to the typical statistical uncertainty on BSAs. 
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2.1 Local fits 



There are at most 12 </)-bins in each Hall B {xB,Q^,t)-hm, and 24 0-bins in each Hall A 
{xB,Q^,t)-hm. ReTi and ImTi are the free parameters of the fits. 

We estimate the systematic errors associated to our if-dominance hypothesis by first 
fitting data setting the subdominant GPDs to 0, then fitting the same data setting the 
subdominant GPDs to their VGG value, and computing the difference. 



2.2 Global fits 

Turning to global fits will help to decrease the statistical uncertainties on the fitted pa- 
rameters. Moreover, we will benefit from the wide kinematic coverage of Hall B data, and 
of the accuracy of Hall A measurements in the same fit. 

Since we are interested only in extracting values of CFFs, we will not try to extrapolate 
outside the kinematic region of the measurements we consider. This allows us to use a 
polynomial parametrisation to perform the fits. The forthcoming difficulty will be the 
evaluation of the systematic uncertainty related to that phenomenological choice. 



2.2.1 A parametrisation of H from the dual model 

The singlet combinations if+ is : 

t, Q') = H{x, t, Q') - H{-x, t, (11) 

This is the quantity which is accessible through DVCS. 

In the framework of the dual model for a spin | target, and assuming if-dominance, 
can be formally expanded according to 



This formal expansion can be resummed as a Gegenbauer polynomial expansion 
[33] and [3l]) : 

oo 

H^{x, e, t, Q') = 2(1 - x') J2 t, Q')Clil,{x) (13) 

n=0 

The coefficients An are defined by : 

where P21 is a Legendre polynomial and Rnp{C,) is a polynomial the degree of which is 
2n + l : 

^nA^) ^ ^) r(n-p + l)r(| + 2p)^ '^V '''2+'' + ^'2+^^'^ j 

(15) 
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with 2F1 the Gauss hypergeometric function. 

The Q^-evolution of Bpi(t, Q"^) is given at leading order in [32] : 



B,,{t,Q-') = B,Ml) ( ) (16) 



where /?o = H - and, for < J9 < N^^^ - 1 ([M]) : 

4 1 



Iv 3 



where we note \l/ the Digamma function and 7^; the Euler-Mascheroni constant. 

We use A = 373 MeV in the MS scheme with 3 flavours of quarks. We obtained this 
value after a running of the strong coupling constant computed at four loops [35] , starting 
from the 2008 world-averaged value of as{Mz) and crossing each quark threshold at (twice) 
its 2008 averaged mass |36j. This evaluation is in good agreement with a recent textbook 
one [37]. The reference scale Qq has been set to Qo = 3 GeV^. 



2.2.2 Iterative fitting procedure 



In practice we truncate Eq. (13) at some maximum value A^max of n and we assume the 
following form for the coefficients Bpi : 

with to a constant and api and hpi the free parameters. Their number is iVmax * (Aniax + 3). 
Due to the truncation at n = A^max the representation of the GPD that we get from a fit 
of {api,bpi) can hardly be trusted outside the domain of the fit. We take it as a smooth 
parametrisatin of the data. 

The selected JLab data consist in 1001 measurements with < ^- We fitted them 
with A'inax = 2, 3, 4 which corresponds to 10, 18 and 28-parameter fits (performed with 
Minuit [38] )• To constrain the polynomial oscillations in ^, we adopt an interative fitting 
procedure. We simplify the problem by first working with bins for which we can neglect the 
i-dependence. Indeed, 40 % of the data with ^ < I satisfy 0.2 GeV^ < -t < 0.4 GeV^. 
We choose to = -0.28 GeV^ in Eq. (18) and extract a first value of the parameters a^/. We 
then add bins with —t between 0.09 and 0.2 GeV^, and between 0.6 and 1.0 GeV^, fitting 
both ttpi and bpi, initialising the fit at the values of ttpi fitted at the previous iteration, 
and bpi at 0. We then add the two last t-bins to the fit, using at each step the previous 
extractions of Qpi and bpi. 



2.2.3 Systematic uncertainties 

The fits with A^max are qualitatively similar. Their x^/d.o.f. are respectively 1.73, 1.61 and 
1.78. The comparison of the values of the CFFs derived from these fits gives an estimate 
of the systematic error on Ti. induced by the truncation. 
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Since we assume if-dominance, we must take into account the systematic error linked 
to the neglect of subdominant GPDs. We proceed as in the case of local fits : we fit the 
data with the subdominant CFFs set to or to their VGG value, and take the difference 
as an estimate of the systematic uncertainty. 



3 Results 

3.1 Extraction of ImH and ReTi 

The Fig. [ijand Fig. |2] display the effect of the truncation of the series in Eq. (13). When ^ 
is small, A^max=2 is not enough to describe the BSAs. On the contrary, when ^ is large, we 
cannot control the fit with A'^max=4. This comes from the fact that statistical errors on BSAs 
are getting larger when ^ grows. The fit with A^max=3 is always good, and close to the local 
fits, which are optimal by construction. We also see that the VGG model overestimates 
the data, which is a known feature [12]. It presumably stems from an overestimation of 
the imaginary parts of CFFs by VGG. At last, global and local fits to helicity-dependent 
and helicity-independent cross sections are all good, and almost indistinguishable. 

We thus choose the fit with A^max=3 as our nominal solution. The systematic uncer- 



tainty on ImTi and ReH linked to the truncation of Eq. ( 13 ) is estimated as the maximum 
of the (absolute values of the) difference between the results for the nominal solution and 
for the 2 other fits. 

The Fig. [Sj Fig. g Fig. [5] and Fig. [6] display our results for IniH and ReTi respec- 
tively. Both local fits and global fit give results with comparable accuracy for ImTi, but 
as expected the results of the global fits are smoother. This is especially true concerning 
ReTi. : in this case the local fits suffer from large fluctuations of ReTi with values which 
fall outside the plot range. However, we could not reliably extract values of the CFFs for 
the larger values of ^ with the global fit. This is reminiscent of the difficulty in controlling 
the oscillating behaviour of the polynomial expansion displayed in Fig. [l| 

The results for local and global fits are almost always compatible, which is a strong 
consistency check. Both rely on the assumptions of twist 2 accuracy and of if-dominance. 
On one hand, local fits suffer from numerical fluctuations (the 2-parameter local fits are 
not constrained enough on some bins) but are almost model-independent. On the other 
hand, global fits are smoother, but suffer from oscillations. That both methods give the 
same results indicates that fluctuations and oscillations are reasonably controlled in the 
bins for which results are displayed. Since, in both cases, the total error bars have the same 
size, we conclude that our estimation of systematic uncertainties due to the truncation of 



the series Eq. ( 13 ) is realistic. 



All fits keep data satisfying ^ < ^- For local fits, changing the maximal value of ^ 
amounts to dropping points. For global fits, the whole results may be changed, but the 
good agreement between the results of both types of fits, and the slow Q^-evolution of the 
extracted CFFs, indicate that this restricted kinematic region is suitable for an analysis 
in the GPD framework. Nevertheless, we observe a sizeable scaling deviation on ImTi. 
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i|) (deg) 360 



Figure 1: The x^-range is divided into 5 bins. The BSAs are displayed at the lowest 
and for the value of t used at the first iteration of the fit. We plot the results of the global 
fits with A^jnax=2, 3, 4, the result of the local fits and the prediction of the VGG code. 
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Figure 2: Helicity-dependent (up) and helicity- independent (down) cross sections ai xb = 
0.36 and = 2.3 GeV^. We plot the results of the global fits with A^max=2, 3, 4, the 
result of the local fits and the prediction of the VGG code. 
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ImH, Local fits 



ImH, Global fit 
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Figure 3: Q^-behaviour of the extracted values of ImTi of local fits (left) and global fit 
(right) on Hall B kinematics. The error bars include both statistics and systematics. 
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ReH, Local fits 



ReH, Global fit 
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Figure 4: Q^-behaviour of the extracted values of ReTi of local fits (left) and global fit 
(right) on Hall B kinematics. The error bars include both statistical and systematic un- 
certainties. 
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Figure 5: Q^-behaviour of the extracted values of ImTi. of local fits (left) and global 
fit (right) on Hall A kinematics. The error bars include both statistical and systematic 
uncertainties. 



13 



Local fits Global fit 









1 




1 

: 1 

, , , 1 ... 1 . 




- 1 

. . . 1 . . 


. 1 





'■4 -0-3,(GeV^) -0.2 -0.4 -0.3 , (g^v^) -0.2 -0.1 



Figure 6: t-behaviour of the extracted values of ImTi. (up) and ReTC (down) of local fits 
(left) and global fit (right) on Hall A kinematics {xb = 0.36 and = 2.3 GeV^). The 
error bars include both statistical and systematic uncertainties. 
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extracted from Hall A data for t = -0.17 GeV^. The choice of the expression of ^ is also 
related to the issue of scaling. Changing the exact expression of ^ to its asymptotic form 
2^'^^ induces differences on the extraction of CFFs, the amplitudes of which depend on 
Xb and t. Deviations are noticeable but results with both expressions of ^ are compatible 
within error bars. This indicates an effect of higher-order power corrections. In view of 
the whole set of data, the conclusion of early Q^-scaling [12| presumably still holds but a 
higher-twist study is needed to make it final. 

Tab. [T]and Tab. |2] summarize our results. Our error bars are dominated by systematic 
effects. Typically we obtain a relative accuracy of 20 to 50 % on ImTC, which is quite 
satisfactory under the assumption of if-dominance and given the statistical accuracy of 
JLab data. On the contrary, ReTi. is still largely undetermined, and is never extracted with 
a precision better than 50 %. But the imaginary parts of CFF might be the quantities of 
prime importance as stressed in : it may be possible to evaluate ReTi (for instance) 
from the knowledge of ImTi on a wide kinematic range through the use of dispersion 
relations. 



Hall B kinematics 


Xb 


t (GeV") 


Q'' (GeV^) 


ImU 


ReU 


0.1342 


-0.1337 


1.1661 


5.90 (0.01) (1.57) (1.01) 


-3.16 (0.04) (1.70) (4.11) 


0.1763 
0.1767 


-0.1346 
-0.1376 


1.3651 
1.5557 


6.91 (0.02) (0.67) (0.72) 
7.69 (0.02) (0.72) (0.72) 


-3.54 (0.03) (3.60) (2.66) 
-4.92 (0.03) (3.80) (2.47) 


0.2350 
0.2377 
0.2460 


-0.1465 
-0.1448 
-0.1442 


1.6453 
1.8895 
2.1641 


5.43 (0.03) (0.34) (0.42) 
5.59 (0.03) (0.35) (0.46) 
5.28 (0.03) (1.08) (0.57) 


-2.75 (0.04) (5.01) (1.41) 
-3.17 (0.04) (4.33) (1.28) 
-3.30 (0.04) (4.18) (1.16) 


0.3205 
0.3215 
0.3213 
0.3211 


-0.1705 
-0.1719 
-0.1743 
-0.1753 


1.9424 
2.2170 
2.5078 
2.7865 


3.33 (0.05) (1.14) (1.64) 
3.46 (0.05) (0.48) (1.74) 
3.55 (0.05) (0.12) (1.81) 
3.64 (0.05) (0.59) (1.89) 


-2.87 (0.04) (7.62) (0.94) 
-2.74 (0.04) (6.66) (0.87) 
-2.66 (0.04) (5.85) (0.82) 
-2.58 (0.04) (5.21) (0.79) 


0.1341 


-0.2840 


1.1678 


5.37 (0.01) (0.60) (1.82) 


-4.50 (0.03) (4.77) (5.48) 


0.1764 
0.1772 


-0.2798 
-0.2819 


1.3680 
1.5653 


6.34 (0.01) (1.52) (1.54) 
6.99 (0.01) (1.97) (1.55) 


-4.21 (0.03) (3.09) (3.39) 
-5.37 (0.03) (3.56) (3.26) 


0.2466 
0.2487 
0.2525 


-0.2842 
-0.2809 
-0.2814 


1.6881 
1.9490 
2.2131 


4.29 (0.02) (0.26) (0.46) 
4.28 (0.02) (0.20) (0.43) 
4.06 (0.02) (0.65) (0.41) 


-2.41 (0.03) (1.68) (0.72) 
-2.70 (0.03) (1.45) (0.62) 
-2.82 (0.03) (1.32) (0.50) 


0.3399 
0.3431 
0.3447 
0.3480 


-0.3062 
-0.3012 
-0.2966 
-0.2942 


1.9930 
2.3060 
2.6372 
2.9706 


3.10 (0.04) (0.66) (0.77) 
3.52 (0.04) (0.36) (0.82) 
3.92 (0.05) (0.49) (0.84) 
4.62 (0.05) (0.36) (0.78) 


-3.48 (0.03) (3.43) (0.06) 
-3.31 (0.03) (2.56) (0.02) 
-3.10 (0.03) (1.84) (0.09) 
-2.82 (0.03) (0.87) (0.15) 


0.1331 


-0.4929 


1.1611 


4.03 (0.01) (0.64) (2.01) 


-4.27 (0.02) (4.81) (6.03) 


see following page . . . 
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Hall B kinematics 


continued from previous page . . . 


Xb 


t (GeV") 


(GeV^) 


ImH 


ReH 


0.1750 
0.1765 


-0.4910 
-0.4909 


1.3580 
1.5611 


4.96 (0.01) (1.85) (2.03) 
5.47 (0.01) (2.37) (2.04) 


-3.98 (0.02) (3.24) (4.14) 
-4.86 (0.02) (3.19) (3.99) 


0.2524 
0.2486 
0.2504 


-0.4875 
-0.4873 
-0.4883 


1.7039 
1.9485 
2.2028 


3.13 (0.02) (0.57) (0.84) 
3.26 (0.02) (0.54) (0.87) 
3.06 (0.02) (0.36) (0.81) 


-1.91 (0.02) (0.47) (0.77) 
-2.21 (0.02) (0.35) (0.83) 
-2.30 (0.02) (0.62) (0.73) 


0.3443 
0.3501 
0.3555 
0.3600 


-0.4964 
-0.4938 
-0.4889 
-0.4854 


2.0062 
2.3282 
2.6851 
3.0455 


1.63 (0.03) (0.34) (0.26) 
1.95 (0.04) (0.37) (0.19) 
2.53 (0.04) (0.48) (0.02) 
3.25 (0.04) (1.06) (0.25) 


-3.81 (0.02) (2.02) (0.41) 
-4.00 (0.02) (1.31) (0.28) 
-4.06 (0.02) (0.43) (0.13) 
-3.99 (0.02) (1.98) (0.03) 


0.1753 


-0.7741 


1.5516 


2.51 (0.01) (0.68) (0.63) 


-1.60 (0.02) (4.81) (1.75) 


0.2493 
0.2476 
0.2494 


-0.7731 
-0.7694 
-0.7689 


1.6847 
1.9394 
2.1990 


1.31 (0.01) (0.46) (0.05) 
1.24 (0.02) (0.54) (0.03) 
0.99 (0.02) (0.71) (0.14) 


-0.52 (0.02) (1.83) (0.25) 
-0.73 (0.02) (1.62) (0.21) 
-0.82 (0.02) (1.59) (0.13) 


0.3516 
0.3597 
0.3607 
0.3582 


-0.7752 
-0.7684 
-0.7623 
-0.7573 


2.0231 
2.3596 
2.7054 
3.0357 


1.39 (0.03) (1.04) (0.68) 
2.09 (0.03) (0.63) (0.71) 
2.20 (0.03) (0.54) (0.56) 
1.83 (0.03) (0.04) (0.39) 


-4.28 (0.02) (1.45) (0.26) 
-4.98 (0.02) (0.97) (0.49) 
-5.18 (0.02) (1.60) (0.66) 
-5.15 (0.02) (2.04) (0.76) 



Table 1: Global fit extraction of ImTi. and ReTi. on Hall B kinematics. The horizontal lines 
gather bins with common {xB,t) and the horizontal double lines bins with common t. The 
errors are in parenthesis, the first one being statistical, the second and third systematic 
(respectively truncation and subdominant CFFs). 



3.2 Discussion 

The Fig. [T] compare our results to a twist 2 model-independent extraction [2H] and an 
extraction with the BMK formalism [12]. Firstly, the use of the GV expressions creates 
important deviations to the latter extraction. Since the extracted combinations of GPDs 
are not the same, we will not make the argument more quantitative. Secondly, we obtained 
results in very good agreement with [28j, but with errors considerably smaller. 

For a given t, and at = 3 GeV^, the Fig. [s] displays the x^-dependence of the 
CFFs we extracted, and compares them to the predictions of the VGG model. To draw 
the authorised region for ImTi. and ReTi. we had to estimate the systematic uncertainty 
induced by the neglect of subdominant CFFs for kinematic configurations for which we 
have no measurements. According to Tab. [l] and Tab. |2| this error weakly depends on 
for fixed t and xb- We averaged it over each {xb, t)-hm, and used it to draw the color band 
of Fig. [Sj Doing so, we enhanced the oscillating behaviour of the contours of the plot by 
adding discontinuities. Since the errors are mostly due to systematic effects, the points in 
the colored domain have the same probability. 
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Global Fit 



M. Guidal 



C. Munoz-Camacho et al. 



7 - 



6 - 



5 - 



4 - 



I I I 



-1 - 



-3 - 



-4 - 



-0-3 t(GeV^)-0-2 



-0.4 -0.3 ,(GeV^)-0.2 -0.4 -0.3 t(GeV^)-0-2 -0.1 



Figure 7: ReTi vs t on Hall A kinematics (0:^=0.36 and = 2.3 GeV^). We compare our 
results (left column) to those of M. Guidal [28] (middle column) and C. Muiioz-Camacho 

et al. [12] (right column). In the latter column 7i+ 2zf^ ^1 + ^T^^ — jm^^^ is plot, and 

not 7i. The error bars include both statistical and systematic uncertainties. 
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t = -0.15GeV^ 



t = -0.30 GeV' 



t = -0.50 GeV^ 



t = -0.80 GeV^ 
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Figure 8: Comparisons of the extracted ImH and ReH and the predictions of the VGG 
model. The error bars include both statistical and systematics contributions. Since our 
errors are dominated by systematics, we can only say that the true value of ImH and ReH 
lie in the colored bands. Moreover, we estimated the systematic errors due to the neglect 
of subdominant GPDs by averaging over the points where measurements are made, hence 
enhancing the oscillating behaviour of the fitting curve (see text for more explanations). 
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Hall A kinematics 


Xb 


t (GeV') 


(GeV^) 


ImH 


ReH 


n '^finn 

u.ouuu 
0.3600 
0.3600 


n 1 7nn 

-U. 1 1 uu 

-0.1700 
-0.1700 


1 .ouuu 
1.9000 
2.3000 


5.87 (0.07) (0.12) (0.57) 
7.05 (0.07) (0.31) (0.45) 


-1.00 ^u.uoj yo.v^j ^u.^i J 
-0.39 (0.03) (1.33) (0.38) 
0.36 (0.03) (0.51) (0.31) 


u.ouuu 

0.3600 
0.3600 


U.ZiOUU 

-0.2300 
-0.2300 


1 . uuuu 

1.9000 
2.3000 


O.UO lU.UUI lU.C/^l lU.OZ,) 

5.83 (0.06) (0.36) (0.74) 
6.67 (0.06) (0.19) (0.64) 


9 77 I'D f\'^\ Q9'l fn 1 9^1 

Z,. 1 1 lU.UO) IZi.c/Z,) lU.XZ,) 

-2.04 (0.03) (1.08) (0.15) 
-1.46 (0.03) (0.32) (0.19) 


0.3600 
0.3600 
0.3600 


-0.2800 
-0.2800 
-0.2800 


1.5000 
1.9000 
2.3000 


4.75 (0.05) (0.99) (0.69) 
5.32 (0.05) (0.43) (0.64) 
5.94 (0.05) (0.20) (0.55) 


-3.45 (0.03) (2.51) (0.16) 
-2.90 (0.03) (0.97) (0.05) 
-2.45 (0.03) (0.20) (0.04) 


0.3600 
0.3600 
0.3600 


-0.3300 
-0.3300 
-0.3300 


1.5000 
1.9000 
2.3000 


4.30 (0.04) (0.95) (0.51) 
4.68 (0.05) (0.52) (0.48) 
5.11 (0.05) (0.15) (0.38) 


-3.86 (0.03) (2.20) (0.40) 
-3.47 (0.03) (0.91) (0.24) 
-3.12 (0.03) (0.15) (0.10) 



Table 2: Global fit extraction of IniH and ReTi on Hall A kinematics. The horizontal lines 
gather bins with common {xB,t) and the horizontal double lines bins with common t. The 
errors are in parenthesis, the first one being statistical, the second and third systematic 
(respectively truncation and subdominant CFFs). 



This being stated, we notice that the extracted values of ImTi. have an x ^-dependence 
similar to that predicted by the VGG model. We again observe that VGG tends to over- 
estimate ImTi.. The extracted CFFs are often but not always compatible with VGG. The 
errors on the extracted ReTi and ImTi. in the hatched zone may be underestimated. We 
already noticed that polynomial oscillations are more severe at large xb, where the data 
are less accurate, as shown in Fig. [T]and Fig. |2j However, from Fig. |3]and Fig. |4]we know 
that the extraction is independent of the type (local or global) of the fit for 0.3 < < 0.4. 
Thus these data should not be rejected but it is necessary to stay cautious when interpreting 
them. 

Conclusions 

Working at leading twist, and assuming iJ-dominance, we extracted ImTi and ReTi. with 
two different methods. The local fits do not benefit from the wide kinematic coverage of 
Hall B data, and are bound to produce rather large statistical uncertainties when fitting 
two independent parameters from an asymmetry measured at (at most) twelve 0-bins. 
On the opposite, the global fit dramatically decreases the systematic uncertainties and 
uses the kinematic coverage to disentangle the contributions of the different CFFs at the 
expense of a fitting Ansatz. The good agreement between the results of both extractions is 
a strong consistency check. It demonstrates that the extraction of ImTi. and ReTi and the 
estimation of the systematics related to this fitting procedure produce relatively model- 
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independent results, at least for xb not too large. This estimation of the systematics will 
be refined but these first results are already encouraging. One immediate advantage is the 
smoothness of the results. This comparison and the weak Q^-dependency of the results 
also validates a posteriori the restriction to kinematic configurations with ^ < ^• 

It is one of the first global fits of measurements from different experiments. In spite 
of the difficulty of controlling polynomial oscillations, the (fitted) behaviour of ImTi. and 
ReTi. versus a;^ is similar to the VGG model prediction. The t-dependence of our global 
fit is also in fair agreeement with that obtained in the first GPD analysis of Hall A data 
in the BMK framework. As expected, we find that the VGG model tends to overestimate 
the physical value of IniH. 

Our results are dominated by systematic uncertainties. Their origin is twofold. Firstly, 
we assumed if-dominance and neglected the contribution of E, E and H. Secondly, the 
Gegenbauer and Legendre expansions induce oscillations in the partial sums of the series. 
We may hope to reduce these systematic uncertainties in the near future using additional 
BSA measurements jH] (unpolarised proton target) and [42j (longitudinaly polarised pro- 
ton target) which will put stronger constraints on the global fits. The extension of our 
procedure to smaller values of is in progress. 

However, our total errors are already of reasonable size, since they are comparable to, 
or smaller than those coming from previous extractions. We typically obtain a 20 to 50 % 
accuracy on ImTi, which is already good regarding our hypothesis of H dominance. 

At last, C. Munoz-Camacho et al. [^2] performed the extraction of a combination 
of GPDs with the BMK formulae. They concluded that Hall A data indicates twist 2 
dominance of DVCS through early Q^-scaling. Using the new GV formalism, we tend to 
come to the same conclusion in this study, but we plan to work out a refined analysis at 
twist 3 in the Wandura-Wilczek approximation to make sure we distinguish between power 
law (higher twist) and logarithmic (evolution) behaviours. 

Note added During the writing of this paper, K. Kumericki and D. Miiller released a 
detailed model-dependent fit of Hermes and JLab DVCS measurements. We will compare 
their results to ours (both present results and the on-going extension to smaller xb) in a 
future study. 
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